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BKL maps and Poincare sections 

Orchidca Maria Lecian^'G 
^Sapienza University of Rome, Physics Department, Piazzale Aldo Mora, 5 - 00185 Rome, Raly 

Cosmological billiards arise as a map of the solution to the Einstein equations, when the most 
general symmetry of the metric tensor is implemented, under the BKL paradigm, for which points 
are spatially decoupled in the asymptotical limit close to the cosmological singularity. 
Cosmological billiards in 4 = 3 + 1 dimensions for the case of pure gravity are analyzed for 
those features, for which the content of Weyl reflections in the BKL maps requires definition of a 
3-dimensional restricted phase space. The role of Poincare sections in these processes is outlined. 
The quantum regime is investigated within this framework: as a result, 1-epoch BKL eras are 
found to be the most probable configuration at which the wavefunctions have to be evaluated; 
furthermore, BKL eras containing n » 1 epochs are shown to be a less probable configuration for 
the wavefunctions. 
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I. INTRODUCTION 

Cosniological billiards are a very powerful description of the solution of the Einstein equations in the case of 
the most general symmetry hypothesized for the metric tensor, as shown in [l[ 0], 3: Hj Hi where the one- 
dimensional BKL map (named after Belinski, Khalatnikov and Lifshitz) was established. In particular, under the 
BKL hypothesis, space gradients are neglected with respect to time derivatives in the asymptotic limit close to the 
cosmological singularity. Here, the solution to the Einstein field equations is mapped to the behavior of a billiard 
ball in a closed domain, which undergoes elastic bounces on the sides of the billiard table, and experiences free-flight 
(according to the geodesies induced by the metric) evolution Q, [3|, 0, 110|- In UM^ UM s-i^d UM^ Hj the geometrical 
interpretation of the statistical analysis was achieved by the definition of 2-variable maps. 
The interest in this description of the cosmological singularity has been renewed by the discovery of new structures 

[n, 0, [13, [a, 0, 0, m. 

In 4 = 3 + 1 dimensions, two kinds of billiard are obtained: the 'big billiard', for which pure gravity is taken 
into account, and for which the original BKL scheme is recovered, and the 'small billiard', which coincides with 1/6 
of the big billiard, and which is defined by the symmetry constraints of the theory |22l |. In this paper, the attention 
will be devoted to the case of 4 = 3 -I- 1 pure gravity. 

In this case, the corresponding cosmological billiard is triangular, and the dynamics of the billiard ball that describes 
the solution to the Einstein equations under the BKL hypothesis is fully understood by means suitable symmctry- 
quoticnting mechanisms, which allow one to schematize the evolution of the billiard ball by means of the so called 
BKL maps. These maps can be defined for the reduced phase space [24I or for the Upper Poincare Half Plane (UPHP) 
|24| . The choice of the map depends on the content of Weyl reflections [23I in the trajectory of the billiard ball 
In fact, the presence of Weyl reflections is not explicitly outlined in the maps for the reduced phase space, and. 



as a result, different subcases have to be considered. These maps, and their content of Weyl reflections, permute 
the role of the Kasner exponents, which represent the projective directions in the asymptotical limit close to the 
cosmological singularity. 

In particular, in this paper, the cases will be analyzed, for which the investigation on the UPHP with a 3-dimensional 
invariant measure has to be followed [23 1231. 



The quantum version of the model can be envisaged according to the choice of boundary conditions that can be 




imposed on the wavefunction, i.e. Dirichlet |28|, [29|, l30|, [31| and Neumann [32ll . [3 311 boundary conditions [3: 
the symmetries that characterize the geoemetrical feataures of the billiard table [33 and of the dynamics [3/ 
[2q , and to the information that has to be looked for [Sj] , [33 . 

One of the greatest achievements of BKL maps is to enclose the continuous dynamics of the billiard chaotic motion 
into a discrete map, which is statistically interpreted as the Poincare return map on the sides of the billiard (i.e. 
the sides of the billiard are considered as a 'special' Poincare section), and on which the geometrical features of the 
billiard table as well as the properties of the dynamics allow to implement several symmetry-quotienting mechanisms. 
From the quantum point of view, on the other hand, the sides of the billiard are imposed boundary conditions, 
according to the quantization procedures. The 3-dimensional energy-shell reduced Liouville measure allows one to 
consider Poincare sections different form the sides of the billiard. The knowledge of the complete dynamics and of the 
symmetry-quotienting mechanisms allows one to extract physical information from these different Poincare sections. 
According to the statistical features of cosmological billiards, the most probable configurations are given by particular 
periodic phenomena on the UPHP. i.e. 1 epoch BKL eras are the most probable ones. Because these periodic con- 
figurations on the UPHP do not correspond to periodic configurations inside the Unit Hyperboloid (UH) [22|, when 
the unquoticnted map is considered, it is necessary to implement the appropriate tools on the UPHP. As a result, the 
configurations corresponding to the classical 1-epoch BKL eras are discovered to be probabilistically favored also in 
the quantum regime. 

The paper is organized as follows. 

In Section |T1 the framework within which cosmological billiards emerge is depicted. 

In Section ini the main features of cosmological billiards are recalled. 

In Section IIIIl the 3-dimensional Liouville measure for special Hamiltonian systems is applied to some problems in 

cosmological billiards, where the content of Weyl refiections on the maps on the UPHP requires the consideration of 

a Poincare section different form the sides of cosmological billiards. 

In Section [IVI the implementation of the quantum regime, and, in particular, of the quantum BKL maps is faced by 

means of the tools acquired in the previous sections. 

In Section |Vl the difference between the geometrical features of the billiard table and the properties of the dynamics 



is stressed within the framework of the analysis developed. Brief concluding remarks end the paper. 



II. BASIC STATEMENTS 

In the case of a general inhomogeneous space-time, close to the cosmological singularity [3^ [401, ^^^ Iwasawa 
decomposition of the metric tensor reads 

d 
gy(x°,x') = ^e-2/3°(-°'-")iV-^(xO^^^)^a^(^o^^z)^ (2.1) 

The relevant degrees of freedom are those expressed by the variables /3, while the other degrees of freedom, encoded 
in the lapse function N and in the shift vector N'^ are considered as 'frozen'. 

As a result, the dynamics of the early universe is described by the chaotic motion of a massless billiard ball in a 
triangular potential well delimited by infinite potential walls, and moving along the geodesies of the metric 

da^ = ^ Gatdp^'d^' - Y, {dn^ - (51 ^/3" ) (2-2) 

a,b=l a \ a / 

According to this description, the bounces against the potential walls are elastic, and the geodesic free-flight evolution 

between each bounce is defined by the Kasner velocities v°-, defined by the Kasner exponents p° as p" = f°/X]a^"' 

which obey the constraints X]aP° = (X]aP°)^ ~ 1- 

In this framework, the cosmological singularity, located at t — > in Eq. p.ip . corresponds to the time variable r — >■ oo, 

where dr = ~dt/ y^. The metric (j2.2p defines a Unit Hyperboloid (UH), inside which the chaotic motion along the 

geodesies takes place. This geodesic motion can be further parameterized according to a radial coordinate p, and the 

remaining degrees of freedom are those encoded by the angular 7*^ variables, which define the surface of the UH, such 

that/3" = p7", 7a7'^ = -l. 

The solution of the Hamiltonian constrainio with respect to the variable A = Inp is equivalent to the projection of the 

geodesic motion onto the surface of the UH. By geometrical transformations, it is possible to project the 7 surface 

on the unit disk, and then on the UPHP. In 4 = 3 -I- 1 dimensions, when pure gravity is considered, the cosmological 

billiard table is triangular, as illustrated in Figure [H where the sides of the billiard table a, b, c are drawn. They are 

defined by the equations a : u ~ ~1, b : u = 0, c : u^ — u + v^ ~ 0, where the variables u and v parameterize the 

UPHP. 

On the UPHP, geodesies are (generalized) half-circles, parametrized by their oriented endpoints. For a 2-dimensional 

system, the complete phase space is 4-dimensional, but, given the symmetries of the problem and the elastic bounces 

of the billiard walls, it is possible to show that the complete phase space is 2-dimensional, given by the oriented 

endpoints of the geodesies. 

The 3 corners of the billiard are due to the directions toward which the 3 anisotropic scale factors reach c». 

For the geometry of the billiard, the succession of trajectories of the billiard ball can be ordered according to the 

corner they take place in, and to the orientation (i.e. according to the oriented endpoints). A trajectory joining any 

two walls is named a (Kasner) epoch, while a set of trajectories taking place inside the same corner of the billiard is 

named a (Kasner) era. As there 3 corners and two possible orientation for each corner (the orientation being defined 

by the first epoch of each era), there are 6 kinds of (Kasner) eras. As geodesies are oriented, epochs are named after 

the corner where they take place, and according to the sides of the billiard they fiy within: in Figure [1] a ba epoch 

is plotted. Accordingly, eras are named after the first epoch they contain: i.e., a set of n epochs taking place within 

the two sides a and b is called an era of the ba type is the first epoch is an oriented trajectory starting from the side 

b and ending on the side a. 

On the UPHP, the geodesies arc circles described by the equation 

v'^ + {u-uof =r^, (2.3) 

where the center uq and the radius r depend on the oriented endpoints u^ and u^ as 

uo = ^^^^ r = 2ii^ (2.4a) 

u = uo + rcos9 v = rsmd (2-4b) 



For a thorough review of the solution of the Hainihonian constraint according to the Misner-Chitre variables, see [41 
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TABLE I. Quotiented maps for the bilhard dynamics. 

These maps act on points defining the epochs and eras of the ba type, for which the variables of the restricted phase space have 

the ranges — oo < u" < —1, u^ > for the epoch map, and —2 < u~ < —1, u+ > for the era map. The maps for the variables 

{u~ ,u^) variables of the restricted phase space are given in [22|. The 1-variable BKL map is recovered by marginalizing the 

variable u~ . 

The maps for the complex variable z = u + iv are given in [26[, by making use of the transformations T{z) = z + 1, S{z) — —^, 

and Ri{z) = —z. The 1-variable BKL map is obtained by evaluating the variable z a.t v = 0, and then by considering the 

oriented endpoint of the corresponding geodesies. 

In Figure[TJ a trajectory of the billiard is sketched, and its oriented endpoints on the u axis of the UPHP are indicated. 
The reduced phase space (m~, u~^) is coordinatized by the oriented endpoints, and consists of several squared 'boxes' 
built of squared 'sub-boxes', which express the length of eras according to the variable u~^. The dynamics of the 'big 
billiard' is therefore described as a 'hopscotch game' on the 'hopscotch court' defined as the reduced phase space 
consisting of the two variables u~ and u"^, where (Kasncr) epochs and (Kasner) eras are defined as points on the 
'hopscotch court'. Even though 6 kinds of eras arc possible, it is possible to chose a preferred orientation, such as the 
ba one, onto which all the kinds of eras can be mapped, according to different symmetry-quotiented mechanisms that 
exchange the role of the Kasner parameters. 



A. BKL maps 



The role of BKL maps is to define different symmetry-quotienting mechanisms for the Kasner Parameters. For the 
anisotropy of space postulated within this framework, indeed, the 3 Kasner parameters are related to the 3 space 
directions. According to the possibility of describing the dynamics on the UPHP by means of the reduced phase space, 
the BKL maps can be stated for the variables (m~, u"*"): in the reduced phase space, the BKL maps act diagonally on 
the two variables. In [22|, the relation between the symmetry-quotiented maps and the complete maps are analyzed 
both for the single- variable version of the transformations, as well as for the double- variable version. In [2g|, the 
features of these maps on the UPHP is examined, and the role of Weyl reflections is outlined. 

All the maps are defined by identifying each era with an n-epoch era of the ba type by identifing the first epochs of 
the era. In the (u^ , u+) reduced phase space, the first epoch of each era is represented as being parametrized by the 
oriented endpoints u^ and u^ in a vertical 'box' delimited by —2 < u~ < —1 and n — 1 < u^ < n. 

The two- variable epoch map and the two- variable era map are recalled in Table |T] both for the restricted phase space 
iu-,u+) and for the UPHP. 

The epoch maps act on the oriented endpoints of the trajectories defining the free-flight evolution of the billiard ball 
and define the change to the new epoch by identifying the two sides of the billiard table which form the corner where 
the epoch takes place. 

The era maps describe the succession of eras by describing the number of epoch in each era via the continued-fraction 
decomposition of the values of the variables {u~ ,u~^) of the first epoch of each era. It identifies the different sides of 
the billiard table when a new era starts. 

All these maps generalize the original BKL map, which is a one-variable map for the auxiliary variable u, identified 
with u"*" in the statistical description. 

According to these maps, for the one- variable map, the continued fraction decomposition of u+ contains the succession 
of BKL eras in the dynamics, and u~ can be interpreted as an initial condition, while, for the two-variable map, the 
variable u~ encodes the 'past' succession BKL eras in the dynamics, according to the action of the map T in Table [H 
The 1-variable BKL epoch map and the BKL 1 variable BKL era-transition map are obtained i) by marginalizing the 
variable m~, such that the BKL u>o map is recovered for the variable u = u~^; ii) by evaluating the geodesies on the 
absolute of the UPHP (i.e. for v = 0), and then by considering the 'rightmost' endpoint , which corresponds to u~^, 
as depicted in Figure [TJ 
Going the other way round, it is possible to 'tile' the UPHP by curvilinear (partially overlapping) domains, where 



epoch of the ba type only are defined, according to the number n, as defined in |26j . 

a. BKL Probabilities Within this picture, probabiHty P„ of an era to contain n epochs is given by the area of 
the pertinent subregion of the reduced phase space according to the invariant measure i i^t-'^u-)^ obtained from the 
reduced form a;(u~, u~^) 

u;iu-,u+) = ^\n^^^^. (2.5) 



and reads 



1 Inlf^. (2.6) 

In 2 n(n+2) *> '' 



From the probabihty P„ in Eq. (|2.6p . P„ is a positive monotonicahy decreasing function of n, such that 1-epoch eras 
are the most probable. Furthermore, as shown on |22| . the dynamics of the billiard exhibits preferred directions, such 
that n-epoch eras are more probable when n is odd. Combining the two information together, one can infer that the 
most probable configuration, also for the unquotiented big billiard, is one given by the periodic succession of 1-epoch 
eras, characterized by the values u^ — — $ and u~^ = (p, where $ is 'Golden Ratio', $ ~ 1, 618, and the 'Small Golden 
Ratio' is defined as + 1 = $. Moreover, because the BKL map is unstable under small perturbations. According 
to the study of the Ljapunov exponents of the problem, the BKL map for the variable u~^ is unstable under small 
perturbations, while the variable u~ is not: a sequence of several 1-epoch eras can be obtained by a small perturbation 
of the periodic configuration but this sequence is a finite one, i.e., after it, an era with n ^ 1 epochs will take place 
sooner or later. 



III. INVARIANT MEASURES AND POINCARE SECTIONS 

In this section, the definition of invariant measures for this model will be revised from j22l | , its connection with the 
definition of a Poincare section for the c Hamiltonian problems will be outlined, and the corresponding definition for 
the chaotic motion inside the unit hyperboloid UH will be proposed. 



A. Energy-shell reduced Liouville measure 

The BKL probability for an n-epoch era to take place obtained by integrating the pertinent subdomain of the 
reduced phase space is originated by the analysis of the invariant measures for Hamiltonian systems. In particular, 
the probability P„ in Eq. (j2.6p holds when a Poincare return map for the dynamics is considered. 
Nonetheless, when the complete geodesic flow is taken into account, i.e. when also the geodesic path is considered, the 
complete phase space is 3-dimensional [221 ■ This way, the (u^,w+) space is sheltered by a 'slab' of varying thickness 
S, with < 5 < a{u~^ ,u~). In other words, the energy-shell reduced Liouville measure fl}^ equals the product of 
the symplectic measure (|2.5p by dS ~ d9 / sinO, with < 9 <tt, such that 

n'--^^ =u{u-,u+)AdS, (3.1) 

and the billiard motion becomes a so-called 'special flow' ^3,, i-e. a combination of uniform motion in the vertical 
S direction, , with discrete 'jumps' in the variables u'^ and u~ . The action of the unquotiented big billiard maps on 
the energy-shell reduced Liouville measure for the unprojected motion inside the UH has been examined in |22| . 

Indeed, the path along the geodesic is parameterized by the angle 9; alternatively, one can choose to parameter- 
ize it according to the value u*, where — 1 < m* < 0, i.e. the projection of the position of the billiard ball on to u axis 
of the UPHP, such that cos 9 = uq — u* . This way, the thickness of the slab defined in [22] reads 

a{u+,u-)=ao{u+,u-) + ^ln l:^;^lZ::\ - (3.2) 



b. Weyl reflections As analyzed in [26[ , the BKL maps defined on the reduced phase space do not encode entirely 
the properties of the transformations, for which Weyl reflections are required. As a direct example, one sees that the 
value cos 6' =\ uq — u* \ implied in the definition of cr(u^, u^) in Eq. p.2|) changes its sign under Weyl reflections on 
the UPHP (as oriented endpoints are considered for the BKL map). 

In the following, the energy-shell reduced Liouville measure will be a fundamental tool to implement the BKL maps 
on the reduced phase space by taking into account the Weyl reflections imposed on the UPHP. As a result, it will 
be possible to consider the symmetry-quotienting mechanisms outlined by the BKL maps for the Kasner coefhcients, 
and to keep track of the Weyl reflections without mismatches. 

Furthermore, the definition of BKL maps as in Subsection III Al allows one develop calculations for the sub-boxes of 
the reduced phase space, where these maps are defined, thus defining clearly the initial conditions and eliminating 
ambiguities. 



B. Poincare sections 

The role of Poincare sections for the description of hyperbolic billiards will be explained in the following paragraphs. 
Poincare sections can be considered when a Poincare return map has to be established, and the physical content of 
the Poincare return maps is better encoded in surfaces different from the sides of the billiard itself. 
Poincare sections will be considered on the UPHP and on the UH. 

c. Poincare Sections on the UPHP It is easily verified that it is not possible to choose any section for the big 
billiard, different from the sides of the billiards, according to which the exact evolution of the big billiard is repro- 
duced. Nevertheless, two strategies are possible: either one i) restricts the ranges of the variables u~ and u+ such 
that they parameterize only one kind of epochs. In the (preferred) case of the ba type, this is achieved by imposing 
—2 < M~ < —1 and < m"*" < oo, or, equivalently, M)on the UPHP, selects only the geodesies which belong to the 
tiling that classifies epochs of the ba type. The details of these domain are recalled in Table HIl 

The vertical S{u*) can be interpreted as a particular Poincare section, through which the geodesies pass, different 
from the sides of the billiard, and through which a Poincare return map can be searched. An example of a Poincare 
section for the return map on the b side of the billiard is given in Figure [TJ 

More in general, one seen that the shape of the Poincare section S is that of a geodesic (i.e. a generalized circle on the 
UPHP defined by u ~ u*) that separates the ca corner of the billiard for the other 2 by exactly cutting the ba corner. 
A Poincare section of this kind is sketched in Figure [T] by the vertical u = const line. This way, one can deduce that 
any geodesies having this behavior, ore even more general ones, can be chosen for these purposes. More in detail, 
any geodesies u = u* with — 1 < u* < 0, or any half-circle intersecting two sides of the billiard can be considered a 
Poincare section for the Hamiltonian flow. Going the other way round, one can interpret the sides of the billiard as 
special Poincare sections (at which the billiard ball bounces) where the return map is considered. 



d. Poincare Sections for the UH The geodesic chaotic motion inside the UH is governed by the metric (|2.2p . 
decomposed along the radial direction p and the angular variable 7, and the (gravitational) walls can be interpreted 
as well as the 'special' (i.e. complete) Poincare sections where the Poincare return map is considered. Projecting the 
dynamics on the 7 surface eliminates one degree of freedom for the model. 

Nevertheless, it is worth remarking that one could think of performing the Poincare returm map for the dynamics 
also inside the /3 hyperboloid by choosing an appropriate Poincare section for this case. Even though the shape of 
the UH can be 'visualized' in a Euclidean space, even though the (gravitational) walls can be represented as planes 
intersecting the UH in a Euclidean space, the geodesic path is given by the metric (|2.2p . This way, and by comparison 
with the considerations for the UPHP of the previous Section, one sees that an appropriate Poincare section inside 
the UH is any (according to the metric (|2.2p ) plane intersecting the light cone where the UH is placed. One easily 
sees that all such planes will have the same role as S{u*). 

In particular, such planes will isolate a corner of the billiard with respect to the other two, and ,in general, it is 
possible that such planes will not bisect any corner of the gravitational walls. 

e. Poincare Sections and BKL maps It is straightforward to verify that the implementation of the BKL symme- 
tries of the dynamics inside the UH, where the Hamiltonian constraint has not been solved yet, play the same role 



of the solution of the Hamiltonian constraint with respect to the variable A (which eliminates one degree of freedom 
by projecting the dynamics on the 7 surface). This would not be true if the BKL syminetries were not implemented, 
because of the 'anisotropic' behavior the complete big billiard exhibits with respect to the existence of a preferred 
direction in the succession of corners. 

As an example, one can choose a generalized (according to the metric (|2.2|) ) plane, different form the gravita- 
tional walls, inside the light cone, of equation /3° = const. The chaotic motion of the point particle inside the UH is 
such that the particle will bounce against the 7 surface. Some trajectories will pass through the section /?" — const 
and some will not. Anyhow, it is possible to identify all these trajectories by identifying the Kasner coefficients (by 
means of suitable Kasner transformations). Because the Kasner coefficients describe the asymptotic (toward the 
singularity) behaviors of the anisotropic scale factors, the Poincare return map of the (now oriented) geodesies on the 
2-dimensional Poincare section (/3° = const in this example) constitute a 2-dimensional statistical description inside 
the UH for the solution of the Einstein equations. 

For this representation of the solution of the Einstein equations, the phase space is 3-dimensional: it consists of 
the two coordinates on the /S"' = const plane, sheltered by a 'slab' which accounts for the direction with which the 
trajectory pass through the generalized plane. 



IV. QUANTUM REGIME 

The quantum implementation of the model is obtained by solving the Hamiltonian constraint, such that a proper 
WDW equation is obtained, and then by considering the quantum properties of the corresponding scheme. The 
quantum regime of this model has been examined both for the UH variables (decomposed as /3° = pj""), and on the 
UPHP (for the variables u and v). 
For the UPHP, the wavefunction ^, defined for the variable z = u + iv oi the UPHP, obeys the eigenvalue equation 

-Alb*(^) = S*(^), (4.1) 

where Alb is the Laplace-Beltrami operator on the UPHP, VP is decomposed along its modes for the variable v, and 
along its Fourier modes for the variable u, such that ^I^ = J^s ^s, where $s is defined as 

$s(m,w) = y^ Cf^y^/^K. _i (27r I ^ I v)exp[2'Kifiu] : (4.2) 

i-i- 

in Eq. (|4.2p . the functions K, are the modified Bessel-X functions of the second kind. As a result, the wavefunction 
^ explores the energy levels s of the eigenvalue equation (j4.1[) . and scans the quantum BKL numbers n through the 
decomposition along the u direction. 

In [26[, the complete wavefunction inside the /3 hyperboloid has been taken onto account; the complete eigenvalue 
equation has been solved by requiring that the wavefunction factorizes into the p and the 7 variables. As a result, 
the wavefunction and its p derivatives have been shown to vanish in the vicinity of the cosmological singularity. 
The boundary conditions for this functions have been examined according to different viewpoints. In particular, 
Neumann and Dirichlet boundary conditions have been chosen in order to stress on different properties. 
Furthermore, in [26[, the symmetry which have to be taken onto account if one wants to interpret the full billiard 
wavefunction as a suitable junction of the wavefunctions of the small billiard have also been considered, according to 
both the symmetry of the billiard table and to the symmetries that govern the dynamics. 



A. Quantum BKL configurations 

Quantum BKL probabilites have been proposed in [26[ , where they have been defined as the integral of the squared 
modulus of the wavefunction VP on the pertinent regions of the UPHP, where the geodesies describing an n-epoch BKL 
era tile the billiard table. Nevertheless, this 'tiling' is not uniform, and these regions overlap. This problem can be 
faced in the reduced phase space (m~, u"''), where the subregions defining n-epoch BKL eras are squared 'boxes' and 
do not overlap. For this purpose, however, one has to consider the full {u~ .,u'^ ,u*) description. Indeed, the Fourier 
decomposition along the u direction of $s is defined for — 1 < u < 0, but detecting the mechanism of quantum BKL 
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TABLE II. The values of the maximum radius rM{n) and of the minimum one rm{n), and of the center uo{n), specified for 
the value n, of the geodesies delimiting the curvilinear domains that tile the UPHP according to the presence of BKL epochs 
of the ba type. Particular kinds of these domains are plotted in Figure [2l 

probabilities implemented by the Poincarc return map on the sides of the billiard can be made difficult by the choice 

of particular boundary conditions for the wavefunction. 

In the following, this inconvenience will be removed by considering the full (m^, u+, u*) description, i.e. by considering 

a suitable Poincare section defined by the generalized geodesies u = u*, with u* not coinciding with any side of the 

billiard. Even though many configurations can be hypothesized, the use of BKL maps recalled in Subsection III Al will 

allow one to implement each problem for the intervals —2 < u~ < —1, u'^ > on the restricted phase space, and 

iJ > 0, -1 < u < on the UPHP. 

The oriented geodesies that define the trajectories of the billiard ball intersect u — u* are parametrized along each 

geodesic path by the variables 



v"^ + {u* - Uof 



(4.3a) 
(4.3b) 



there are only 2 equations for the set of three variables (u^,u^,u*), and this is required for the definition of the 
position of the billiard ball along a given geodesies if the variable u has to be decomposed according to its Fourier 
modes in the definition of Maass wavefunctions. 
As a result, the eigenfunctions for this model will read 



$ 



:{u,v) = (r^ - {u* - UQ)'^)''-/'''^^Cf^exp[7Tifj.u*]IC^_i{2Tr \ fi \ {\Jr'^ - (u* - uq)^)) 



(4.4) 



As outlined in [26| , on the UPHP, the billiard table is tiled according to the number n of epochs contained in BKL 
eras, where the tiles are delimited by the geodesies Q^^ in) and Q™'(n), whose radius and centers are recalled in Table 
im Nevertheless, this tiling is not uniform, and the corresponding subregions partially overlap, as shown in Figure [H 
For this, given two different geodesies parameterizing an n-epoch era and an n + 1-epoch era, respectively, it is 
possible to find a configuration for which they intersect at a given — 1 < w* < and a given v. The difference in the 
number of BKL epochs of the era they represent is given by the direction they pass the Poincare section m = u*, i.e. 
by the angle Q, in Eq. (|3.2p . This information is not entirely encoded in the oriented endpoints of the geodesies u~ 
and u+. 

When the quantum wavefunction is considered, two different wavefunctions $s,ti and $s_„+i acquire the same value, 
and so do their derivatives with respect to v. the only difference is given by the quantum BKL numbers n and n + 1, 
and this information is encoded in the derivatives with respect to w"*" and u^ . 



To investigate the physical information enclosed in the wavefunction $s, it is convenient to define the object 



A^^, (u ,u+,u , u+ ) as the ratio between two different derivatives of the wavefunctions, evaluated at the same 
point u = u^ and v = v(u*) (through which an infinite number of geodesies pass, with different endpoints), which 
encode different successions of BKL epochs and eras, as 



^nt(.U-,U+,U-',U+') 









(4.5) 



when evaluated at two different endpoints u and u , it encodes the different successions of BKL epochs and eras, 

as for their continued-fraction decomposition and for the action of the map T. 

The interpretation of the one dimensional BKL map, i.e. that where only the value u"*" is taken into account, consists 

in considering the values of u~ as 'initial conditions'. 

For the case m+ = m+ + 1, for large values of n, one obtains 



A++ ^ 1 



(4.6) 



where m is the value of epochs contained in the previous BKL era, according to the definition of the map T. Eq. 
(|4.6p has the same asymptotical (for large n) series expansion of the ratio between the two BKL probabilities Eq. 
(|2.6p evaluated at n and at n + 1, i.e. Pn/Pn+i — 1 + 1/n^. This allows one to infer that the ratio between the two 
derivatives, evaluated at the same point u* , and with the same boundary and initial conditions, states the asymptotical 
probability to find an n-epoch BKL era according to the same classical probability law (|2.6p . for which BKL eras 
containing a number n >> 1 of BKL epochs are probabilistically disfavored, and, in particular, n-epoch eras are more 
probable than n + 1-epoch eras. 
For any value of n and n', one has 

A _ «+--»* _ n-l+x+-u* . /^ -N 

n,n+l u+'—u' ^ n' — l+x'. —11." ' \^- ' ) 



this ratio is extremized for small values of n, i.e. for n = 1, at a fixed n' ^ n. In Eq. (|4.7p . x+ is the fractional part 
of u^ . 

By the action of the map T on u~ , one then learns that the most probable configuration, also form the quan- 
tum point of view, is one for which u^ is close to the value of $. 

As a result, it is possible to conclude that, also from the quantum point of view, the most probable configura- 
tion at which the wavefunction $s has to be evaluated is (at least) close to that corresponding to the periodic 
configuration of 1-epoch eras. 

The point in having introduced u* is avoiding to establish any fixed boundary condition (such as Neumann or 
Dirichlet boundary conditions) on the wavefunction. In fact, it is in principle always possible to define a Poincare 
section on the billiard table, where the wavefunctions do not vanish identically. 



B. Quantum BKL probabilities 

Collecting all the ingredients together, it is possible to express the quantum BKL probabilities by means of a 
suitable integral on the pertinent regions of the reduced phase space. It is interesting to recall that it is not possible 
to define a Poincare section a priori on the restricted phase space {u~ , u+), because the energy-shell reduced Liouville 
measure p.ip is composed by the two elements ui{u~ ,u^) and dS as 'complementary' (or orthogonal). 
The formal expression for the quantum BKL probability P„.s is then given by 



du 



n-l 



du^ £-eT^^^\^s{u~.u\u*)\\ (4.8) 



where the integration is performed according to the classical energy-shell reduced Liouville measure (|3.ip recalled 

in Section IlIII It is important to remark that here the BKL era map is being used, such that the orientation of the 

angle is always the same, and no sign ambiguities arise. Integrating the geodesic path dS{u*) along a complete 

BKL (oriented) trajectory consists in considering the two values u* = —1 and u* = 0, for which cos6'i =| mq + 1 | and 

cos 02 =1 "0 I (as one can easily verify from Figure [H). 

With respect to the formal result obtained in [2g| , where the integration was performed on the UPHP on overlapping 

domains, this way the problem is avoided by following the geodesic path inside the billiard and then by integrating 

the remaining functions on the reduced phase space. 

Eq. (|4.8p defines the quantum BKL probabilities for a BKL era to consists of n BKL epochs form the quantum point 

of view. 

It is possible to express also the probability for an epoch to be the n'-th epoch of an iV > n'-cpoch BKL era, 
as 

-2-n' JN-n' J 02 

where the quotientcd 7f,a map has been used for the evaluation of the integration boundaries, and the notation of [26[ 
has been kept for the labellings of epochs. 
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V. OUTLOOK 

In this analysis, some features of the cosmological billiards describing pure gravity in 4 = 3 + 1 dimensions have 
been investigated. Particular attention has been devoted to the difference between the geometrical properties of the 
billiard domains and the dynamical properties of the billiard chaotic evolution. In fact, when symmetry quotienting 
mechanisms have to be implemented, such that BKL maps can established, some differences between the two proce- 
dures arise, and they can be interpreted, at least in the most straightforward analysis, according to the content of 
Weyl reflections in these maps: even though it is, in principle, possible to identify the sides of the billiard without 
Weyl reflections, the information needed to consistently encode the dynamics requires this kind of reflections. 
These only apparent 'mismatch' is due to the features of the Einstein field equations, where the 3 Kasner directions 
are considered as an unordered set, while the symmetry-quotienting mechanisms are based on the orientation of the 
geodesies. 

A possible strategy for reconciling these different viewpoints is to consider the 3-dimensional energy-shell reduced 
Liouville measure fl)^ in p.ip . by which the content of Weyl reflections in the UPHP can be taken care of. 

From the quantum point of view, the difference between imposing boundary conditions for the wavcfunction ac- 
cording to the features of the billiard table and to those of the BKL maps for the dynamics has been investigated 
in [], and the deflnition of quantum BKL probabilities has been formally proposed. This definition can be made 
more effective, such that more physical information can be extracted from it, if the wavcfunction is expressed in the 
{u~ ,u~^) variables of the restricted phase space by introducing a suitable Poincare section on the UPHP, and then by 

(3) 

formally integrate it by means of the invariant measure i^)^ ■ By this procedure, it has been possible to investigate 
what quantum BKL numbers are probabilistically favored. Interpreting the available information according to the 
continued- fraction decomposition of the (u^, m^) variables, is is possible to infer that 1 epoch BKL eras are the most 
probable also in the quantum regime. 

Because Poincare sections of the big billiard (i.e. on the UPHP) can be verified to contain by construction dif- 
ferent kinds of eras (i.e. not only of the ba type), the action of any symmetry-quotienting mechanism requires further 
specifications, as it has been necessary to impose throughout the present analysis. 

From the quantum point of view, these procedures allow one to gain physical insight about the BKL configurations 
describing the succession of epochs and eras without imposing any specific boundary conditions. Indeed, the statistical 
meaning of BKL maps can be researched for in the Poincare return map of the bounces of the billiard ball on the 
sides of the billiard table, on which boundary conditions arc usually imposed boundary conditions on. 



VI. CONCLUDING REMARKS 

In this paper, some particular aspects of the cosmological billiards describing pure gravity in 4 = 3 + 1 space time 
dimensions have been investigated, and some light has been shed both at the classical level and as far as the quantum 
regime is concerned. In particular, the phenomena have been analyzed, for which the content of Weyl reflections on 
the maps on the UPHP requires the introduction of a 3-dimensional invariant measure, when physical objects have 
to be defined by means of integration of the corresponding invariants on suitable regions of the phase space. 
From the classical point of view, the implementation of BKL maps on for the motion inside the UH has been consid- 
ered as far as the action of BKL in exchanging the role of the Kasner directions is concerned. 

Form the quantum point of view, the definition of quantum BKL probabilities for the BKL quantum numbers con- 
tained in the wavefunctions solving the eigenvalue problem on the UPHP has been refined by the formal result (|4.8p . 
In this case, the squared modulus of the wavcfunction has to be integrated by means of the 3-dimensional energy-shell 
restricted Liouville measure. 



Furthermore, by analyzing the behavior of the wavcfunction via the tools refined in |22| and |26|, by the comparison 
of the derivatives of the wavefunctions with respect to the value of the oriented cndpoints which encode the the 
succession of BKL epochs in the dynamics, the most probable configuration form the quantum point of view has 
been defined to the still that of a sequence of 1-epoch BKL eras. In particular, this can be achieved by the periodic 
configuration defined by the values u~ = — $ and u~^ = 0, or by a small perturbation of it. 

More in detail, it is possible to infer this information by interpreting the values of the oriented cndpoints u+ and 
u~ as encoding the number of BKL epochs in the future evolution and in the past one, according to the continued 
fraction decomposition of these parameters, respectively. 
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The paper has been organized as follows. 

In Section |lll the features of cosmofogical bihiards and their description by means of suitable symmetry-quotienting 
maps (the BKL maps) have been recalled. 

In Section IIIIl the physical problems for which a 2-dimensional scenario requires a 3-dimensional phase space have 
been outlined. The use of Poincare sections in the definition of Poincare return maps has been adopted for the 
description of particular features of cosmological billiards. 

In Section llVl the analysis of quantum BKL numbers has been refined according to these tools, and the formal def- 
inition of quantum BKL probabilities has been rested according to the 3-dimensional energy-shell reduced Liouville 

(3) 

measure ^)^ ■ As a result, the most probable configuration on which the wavefunction has to be calculated is one 
close given by the periodic succession of 1-epoch eras, i.e. for n ~ I also form the quantum point of view. 
In Section |Vl the role of the invariant measure 57^ in reconciling geometrical features of the billiard tables with that 
of the dynamics, which is encoded by a suitable composition of Weyl reflections, has been clarified. 
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FIG. 1. An oriented geodesies on the UPHP is plotted in green (light gray): the billiard trajectories are parameterized according 
to the oriented endpoints u^ and it"; the radius of the corresponding half-circle is plotted in blue (dark gray), and the center 
uo is indicated. A Poincare section m = it* is plotted in gray. 
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